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1 INTRODUCTION 



General Relativity (GR) is a beautiful scheme for describing the gravitational 
field. This theory is believed to apply to all forms of interactions, especially 
between large scale gravitational structures. It has been proven that black 
holes exist on the basis of study of the effects they exert on their surround- 
ings. They greatly affect the surrounding plasma medium with their enor- 
mous gravitational fields. Since all compact objects have strong gravitational 
fields near their surfaces [I], it is important to study the general relativistic 
effects on physical processes, like electromagnetic processes taking place in 
their vicinity. Magnetohydrodynamics (MHD) with the effects of gravity is 
called general relativistic magnetohydrodynamics (GRMHD). The GRMHD 
equations help us to study stationary configurations and dynamic evolution 
of a conducting fluid in a magnetosphere. They include Maxwell's equations, 
Ohm's law and mass, momentum, and energy conservation equations. These 
equations are required to investigate various aspects of the interaction of 
relativistic gravity with a plasma's magnetic field. 

The 3+1 formalism is well-suited to carry non-relativistic intuition of 
physicists about plasmas, hydrodynamics, and stellar dynamics into the 
arena of black holes and general relativity. This formalism (also called the 
ADM formalism) was originally developed by Arnowitt et al. |2j. It was 
motivated by several startling results proved in the 1970s by using a black- 
hole viewpoint [3-9]. Thorne and Macdonald [10] and Thorne et al. [UJ 
extended this formulation to electromagnetic fields in black-hole theory. The 
wave propagation theory in Friedmann universe was investigated in a 3+1 
formalism by Holcomb and Tajima [T2] , Holcomb [32], and Dettmann et al. 
[13] . Khanna [T5] derived the MHD equations describing the two-component 
plasma theory of a Kerr black hole in this split. Komissarov [16] discussed the 
Blandford-Znajek monopole solution by using the 3+1 formalism in black- 
hole electrodynamics. Zhang [17,18] formulated the black-hole theory for 
stationary symmetric GRMHD by using the 3+1 formalism with its applica- 
tions on a Kerr geometry. 

The formalism for gravitational perturbations away from the Schwarzschild 
background has been developed by Regge and Wheeler [19] . It was extended 
by Zerilli [20], who showed that the perturbations corresponding to changes 
in the mass, the angular momentum, and the charge of the Schwarzschild 
black hole are well-behaved. The decay of non-well-behaved perturbations 
has been investigated by Price [21]. The quasi-static electric problem was 



2 



solved by Hanni and Ruffini [5], who proved that the lines of force diverge 
at the horizon for an observer at infinity. Wald [22] derived the solution 
for the electromagnetic field occurring when a stationary axisymmetric black 
hole is placed in a uniform magnetic field aligned along the symmetry axis 
of the black hole. A linearized treatment of plasma waves for a special rela- 
tivistic formulation of the Schwarzschild black hole was developed by Sakai 
and Kawata [23]- Buzzi et al. [24] extended this treatment to waves propa- 
gating in the radial direction in a general-relativistic two-component plasma 
by using the 3+1 ADM formalism. They investigated wave modes for one- 
dimensional radial propagation of transverse and longitudinal waves close to 
the Schwarzschild horizon. 

The aim of this paper is to investigate the properties of isothermal plasma 
waves by taking the Schwarzschild spacetime in a planar analogue. The 
scheme of the paper is as follows: The next section provides comprehensive 
details about the results in the 3+1 formalism. In Section III, we obtain 
the GRMHD equations for a perfect fluid with perfect MHD flow conditions. 
Section IV describes the perturbations and relative assumptions that will be 
used in the ensuing sections. In Section V, we apply the perturbations and 
Fourier- analysis methods to the GRMHD equations. Section VI investigates 
the dispersion relations obtained in graphical forms. The last section contains 
a summary and discussion of the results. 



2 3+1 SPACETIME MODELING 

In the 3 + 1 formalism, the line element of the spacetime can be written as 

m 

ds 2 = -a 2 dt 2 + 7ii {dx l + /3 l dt)(dx J + (3 3 dt), (2.1) 

where a is the lapse function, f3 l are components of the shift vector, and 7^ 
are components of the spatial metric. All these quantities are functions of 
coordinates t and x l . A natural observer, associated with this spacetime and 
called the fiducial observer (FIDO), has a four- velocity n perpendicular to 
the hypersurfaces of constant time t and is given by 

n = "Oar ~P 7TT ' 
a at ax 1 

Notice that we are using geometrized units so that G = 1 = c. In four- 
dimensional spacetime, vectors and tensors will be denoted by boldface italic 
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letters. A letter having a dyad over it will show a three-dimensional tensor. 
All vector analysis notations, such as gradient, curl and vector product, will 
be those of the three-dimensional absolute space with the three metric 7 . 
The Latin letters k, ... represent indices in absolute space and run from 
1 to 3. 

The perfect MHD flow assumption [T7] is 

E + VxB = 0, (2.2) 

where V, E, and B are the velocity and the electric and magnetic fields of 
the fluid as measured by the FIDO. This shows that there can be no electric 
field in the fluid's rest-frame. 

Applying this condition, Faraday's law in the 3+1 formalism becomes 

— + — (B.V)/9 + 9~B = —V x (aV x B), (2.3) 
ar a a 

where the FIDO's four-velocity expansion rate and FIDO's measured rate of 
change of any three-dimensional vector in absolute space (i.e., orthogonal to 
n) can be expressed, respectively, as 

0=i(f-V./J), jUi(|-/3.V). 
a 2g ar a at 

For perfect MHD flow, the equations of evolution of the magnetic field give 

DB 



+ -(B.V)(/3 -aV) + <9 + -V.(aV) } B = 0, (2.4) 



Dt a |a 

where 

is the time derivative moving with the fluid. The local conservation law of 
the rest mass according to the FIDO is 

^WV.^{g + V.(aV-^0, (,5> 

with po being the rest-mass density. The law of force balance measured by 
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the FIDO is given by 

2 B 2 \ 4 1 1 DV j 2 Da 

Pol V + J^J Tii + Po7 A*W - ^BiBjj — + Po7 ^— 

-PV + -L(V x B) 4 V.(V x B) - J-,(«B)5 + -^(a*)^ 

1 -[B x {V x (V x (qV x B) - (B.V)/3) + (V x B).V/?}]<, (2.6) 



Ana 



where p is the pressure. The FIDO's measured local energy conservation law 
[T7] is given by 

/ + fle + — lr'-'U'r.'M = ^V.(« 2 S) + -V(3 : W + E.j. (2.7) 

ar 2a a z a 

When we substitute the following values for a perfect fluid, 

e = ( m -p(l-V 2 )) 7 2 , 
S = PPo7 2v , 
W = /ipo7 2 V® V + pY, 

with /i being the specific enthalpy, and perfect MHD flow assumption, Eq. 
(12.71) becomes 

Po7 2 ^ + Pl 2 ^ + 2p p 7 4 V.^ --£+ ( m 7 2 - P)9 + 2p 0/ u 7 2 V.a 
+p 0/ u 7 2 (V.V) - i{ Po p 7 2 V.(V.V) / 3 + p(V./3)} + ^{p P7 2 ^^ 

+Pl ij }£tH j + t^-[(V x B).(V x (aB)) + 6»a(V x B) 
Ana 

+a(V x B).— (V x B) + (V x B).(V x B.V)/3.(V xB)]=fl (2.8) 
ar 

with acceleration a = — . Equations (I2.3p - (I2.6I) and (I2.8P are the perfect 
GRMHD equations for a plasma in the vicinity of a black hole. 
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3 GRMHD EQUATIONS FOR ISOTHER- 
MAL PLASMA IN SCHWARZSCHILD PLA- 
NAR ANALOGUE 



This section will give a detailed structure of the chosen spacetime with the 
respective form of the GRMHD equations. 

3.1 Description of the Planar Schwarzschild Spacetime 
and the Equation of State 

The planar approximation of the Schwarzschild line element [18] is 

ds 2 = -a 2 (z)dt 2 + dx 2 + dy 2 + dz 2 , (3.1) 

which is an analogue of the Schwarzschild spacetime with z, x and y as radial 
r, axial 4>, and poloidal 9 directions, respectively. The lapse function vanishes 
at the horizon, which can be placed at z = 0, and increases monotonically 
to unity as z increases from to oo. The relationship between the universal 
time t and the FIDO's time r can be expressed by the time lapse a = 
The Schwarzschild black hole is non-rotating; hence, there is no shift of 
coordinates. In this spacetime, the FIDO has a four-velocity equivalent to 

a at 

In the hydrodynamic treatment, the plasma is represented as a perfect 
fluid. In a planar analogue of the Schwarzschild spacetime, we assume that 
the system of perfect GRMHD equations is enclosed by the isothermal state 
of a plasma. This state can be expressed by the following equation [17J: 

_ P + V _ cons i an i^ (3.2) 
Po 

which shows that there is no exchange of the heat between the plasma and 
the magnetic field of fluid. 

3.2 Respective GRMHD Equations 

For the Schwarzschild planar analogue given in Eq. (13.ip . the GRMHD 
equations Eqs. fl2.3l) - fl2.6p and (12. 8p . in the 3+1 formalism for the isothermal 
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state of a plasma take the form 

_ = V x («V x B), (3.3) 

— + (aV.V)B - (B.V)(aV) + BV.(aV) = 0, (3.4) 
^ + ( £ ,V.V)(, + P ) + (p + P hYf 

+(p + p)7 2 V.(aV.V)V + (p + p)V.(aV) = 0, (3.5) 

{ ^ + ^ + 9 ^ + + ~ h* B > } K + V - V ) ^ 

- " ^ = -(P + P)7 2 a, + ^V.(V x B)(V x B), 
(MB) " (nC ' ) -' L> ' ' [Bx{Vx(Vx(aVxB))}] i -p, i) (3.6) 



8na 2 Aira Ana 

7 2 (^ + V.V)(, + p) + 2(p + ^V.(i| + V.V)V 

+2(p + p) 7 2 V.a - -% + (p + p) 7 2 V.V 
a at 

+-^-[(V x B).(V x (aB)) + a(V x B).-^-(V x B)] = 0. (3.7) 
Ana a at 

Equations fl3.3l) - fl3.7p are the perfect GRMHD equations for an isothermal 
plasma in the vicinity of a Schwarzschild black hole. In the rest of the paper, 
we shall analyze these equations by using perturbations and Fourier analysis 
procedures. 



4 RELATIVE ASSUMPTIONS 

In this background, the perturbed flow of the fluid is considered to be along 
the z-axis. Thus the FIDO's measured four-velocity and magnetic field are 
along the z-axis and can be expressed as 

V = u(z)e z , B = Be z , 

where B is a constant quantity. For this velocity, the Lorentz factor 7 takes 
the form 7 — ! 



VT- 
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The perturbed flow in the magnetosphere shall be characterized by its 
fluid density p, pressure p, velocity V, and magnetic field B (as measured by 
the FIDO). The first-order perturbations in the above-mentioned quantities 
are denoted by 5p, 8p, 5V, and <5B. Consequently, the perturbed variables 
will take the following form: 

p = p° + Sp = p° + pp, p = p° + Sp = p° + pp, 
B = B° + = B° + Bb, V = V° + 5V = V° + v, (4.1) 

where p°, p°, B° and V° are unperturbed quantities. Due to gravitation, 
the waves can propagate in the ^-direction with respect to time t. Hence, 
the perturbed quantities depend on z and t. 

The following notations for the perturbed quantities will be used: 

p = — = p(t,z), p 
P 

v = 5V — v z (t, z)e x , b 

We also assume that the perturbations have harmonic space and time de- 
pendences, i.e., 

p, p, vandb- e -*(«*-to). 
Thus, the perturbed variables can be expressed as 



p(t,z) 


= c\e 


-i(ujt- 


-kz) 

; 


p(t,z) 


= c 2 e" 


-i(u)t- 


-kz) 

; 


v z (t,z) 


= c 3 e" 


-l(il)t- 


-kz) 


b z (t,z) 


= c^e' 


-t(uit- 


-kz) 



— =P(t,z), 

p 

5B 
~B 



b z (t, z)e z 
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5 PERTURBATION EQUATIONS 

If the perturbations given by Eq. (14. ip are introduced, the GRMHD equa- 
tions, Eqs. (I3.3p - fl3.7l) . become 

= V x (ov x B) + V x (qV x SB), (5.1) 
V.(5B) = 0, (5.2) 

d ^ P dt 6 ^ + ^ aV ' V ^ 6p + ^ + ( p + P ^ Y ^ ~ + ^ V ' V ln u 
+a(p + p)(V.v) + (5p + 5p)V.(aV) + (6p + 5p) 7 2 V.(aV.V)V 
+2(p + p) 7 2 (V.v)(aV.V) ln 7 + (p + p) 7 2 (aV.VV).v 
+ (p + p) 7 2 V.(aV.V)v = 0, (5.3) 
R 2 1 1 8vi 

{((p+p)7 2 + ^ + (P+Ph'm - ^B iBj }-— 

+i_[ B x {V x + (p+ph'v^Vi + {p + ptfViVteV'V 

--^-{(at&Bi)^ - (a5Bj)£Bi = - 7 2 {(5p + 8p) + 2(p + p) 7 2 (V.v)H 

~(Sp),i + ^{Milj - (aB^^SBi - (p + p) 7 4 (^' + v 3 Vl )V k)3 V k 
47ra 

- 7 2 {(5p + 5p)V j + 2(p + p) 7 2 (V.v)\/ J + (p + pK' 

- 7 Vi{(5p + <5p) V J + 4(p + p) 7 2 (V.v)V J + (p + pV MvfcV* (5.4) 

^1 ^p + *p) + 2 4y * _ 2(p + ph4(v . v)(v . v)ln , 

a at a at 

+6(p + p) 7 6 (V.v){V.(V.V)V} + (p + p) 7 4 V.(v.V)V 

+2(p + p) 7 4 V.(V.V)v + 2(5p + 5p) 7 2 V.a + 2(p + p) 7 4 (V.v)(V.a) 

2/V7 \ , / , \ 2 „ 2/ T r rrVJt i X \ 1 5 (<^) 



+(p + p) 7 2 (V.v) + (p + p) 7 Va + 7 2 (V.V)(5p + £p) 

a at 

-(p + p) 7 2 (v.V) In w + 2(5p + 5p) 7 4 V.(V.V)V + (8p + <5p) 7 2 (V.V) 

+2(p + p) 7 4 (V.v)(V.V) + 2(p + p) 7 4 v.(V.V)V + J_[(V x B). 

47ra 

(V x (a5B)) + (v x B).(V x (aB)) + (V x <5B).(V x (aB)) 

dSB <9v 
+ (V x B).(V x _) + (V x B).(^ x B)] = 0. (5.5) 
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The component forms of Eqs. (I5.1l) - (l5.5p are given as 

19h 0, (5.6) 



a dt 

b z , z = 0, (5.7) 

d p dip du % 

p ~dt + ua ( p ?' z + pp ~' z ^ + p ~dt + ^ u ^ p + ^ ~dt 
+a(p + p)(l + j 2 u 2 )v z>z — (p — p){(aup)' + a^ 2 u 2 pu} 

= a{p + p){l - 2 7 V)(1 + 7 V)% Z , (5.8) 

(P + P)7 2 (l + l 2 u 2 ) l~~^ + UVz > x j +PP,z+PP ! z 

= -1 2 (P + P)W(1 +7VXI + 47V) + 2u^ 2 a z )}v z 

-(pp + pp)i 2 {a z + u(l + ^ 2 u 2 )u'}, (5.9) 

m ~7^P + p ^ -7nP+-(P + Pn u ^r + {p + P)l (1 + 27 m )v z>z 

a at a at a at 

+ (pp + PP)l 2 {2ua z + (1 + 2^ 2 u 2 )u'} + i 2 u(p^ z p + p tZ p + p iZ p + p tZ p) 

+ (P + P)7 2 l^l 2 uu + a z )(l + 2 7 V) - - I v z = 0. (5.10) 

The conservation law of rest-mass p2] in a three-dimensional hypersurface 
for the isothermal state of a plasma, 

a(p + p)ju = constant, 

is used to obtain Eqs. (I5.8P and (I5.10p . The same law will be applied for 
further simplifications. 

Using Eq. (14. 2[) in Eqs. fl5.6p - fl5.10p . we obtain the Fourier- analyzed form 



10 



of the above equations as 

LOJ 

c 4 = 0, 5.11 

a 

tkci = 0, (5.12) 
c\{— piuo + ikpau — (nap)' — au 2r y 2 pu'} + C2-{ — pioo + ikpau 
+ (uap)' + au 2r y 2 pu'} + c^{p + p) {a(l + ^ 2 u 2 )ik 



u' 



= 0, (5.13) 
cip^ 2 {a z + uu'(l + 7 2 w 2 )} + c 2 {p-y 2 {a z + uu'(l + 7 2 u 2 )} + ikp + p'} 



-ail - 27V)(1 + 7 u ) 2^7^ 

it 



+c 3 (p + ph 2 

+2ua z }} = 0, 



(1 + 7 V) 



-LUJ 



a 



+ Luk) + {u'(l + 7 V) (1 + 4 7 V) 



(5.14) 



PI 



2 / ~ iuJ 
a 



+ uku + 7 up + p27 wa 2 + p(l + 27 m )7 ii > ci 



-to; 



a 



p(7 — 1) + ik'-f up + 7 -up' + 27 wpa 2 + p(l + 27 u )7 -u' > c 2 



+ (P + P)7 2 



-2 to; 



a 



7% + + 2 7 V) + 37W(1 + 2 7 V) 



+a 2 (l + 2 7 V) - - } c 3 = 0. 



(5.15) 



Equations (15. lip and (15. 12[) yield that C4 = 0; i.e., the magnetic field is not 
affected by the black-hole gravity and time. It is mentioned here that Eqs. 
( I5.13p - (|5.15p also give a non-magnetized plasma. 



6 NUMERICAL SOLUTIONS 

We consider that the magnetosphere is filled with a stiff fluid for which p = 
constant = p and let a = tanh(z). Using these values, the mass conservation 
law in 3-dimensions gives u = —, — = . 

yi+tanh 2 (z) 

When we use these values, the determinant of the coefficients of constants 
ci, C2, and C3 in Eqs. (I5.13p - (I5.15I) give a complex dispersion equation of the 
type 

Ai{z)k 2 + A 2 (z, u)k + A 3 (z, to) + L{A 4 (z)k 3 

+A 5 (z, to)k 2 + Ae(z, to) + A 7 (z, to)k} = 0, (6.1) 
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Surface Plot of Propagation Factor Surface Plot of Attenuation Factor 




Figure 1: Waves propagate more near the event horizon. They damp as they 
move away from the event horizon. Anomalous dispersion is found in that 
region. 

which, on solving, gives three complex values of k = k\ + ik 2 - The sinusoidal 
expressions then take the form e -<-M- fc i z - tfc 22) _ e -t(u>t-kiz)-k 2 z _ j g ODV j ous 

that Re(k) is the propagation factor and Im(k) is the attenuation factor for 
a time harmonic plane wave with fixed angular frequency u in a dispersive 
material. Using the values of k, the phase velocity (v p = jjr) and the group 
velocity (v g — [25J of the waves can be calculated. The three values of 
the dispersion relation are shown in Figs. 1, 2 and 3. 

Fig. 1 shows that the propagation factor is infinite at the horizon, which 
indicates that no wave exists there. This has greater value near the event 
horizon and decreases as we go away from the horizon. It also increases with 
increasing angular frequency of the wave. The attenuation factor increases, 
which shows that the waves damp with increasing angular frequency and 
value of z. The phase velocity is greater than the group velocity, which 
shows that dispersion is normal in the region. The waves with negligible 
angular frequency have phase velocities less than the group velocities and 
hence, the medium shows anomalous dispersion. 
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Figure 2: Waves propagate rapidly as they move near the event horizon. The 
dispersion of waves is anomalous in that region. 

The propagation factor in Fig. 2 shows the same behavior as that in 
Fig. 1. In the region, the attenuation factor takes very small values, which 
are positive and negative randomly. The phase and the group velocities are 
the same in the region, except for zero angular frequency waves. For those 
waves, the phase velocity is less than the group velocity, indicating anomalous 
dispersion. 

In Fig. 3, the propagation factor is greater near the event horizon and 
decreases as we go away from the horizon, but becomes infinite at the event 
horizon. The attenuation factor increases when the angular frequency in- 
creases and decreases when the value of z increases. The waves grow as they 
go far from the event horizon of the black hole. Damping occurs when the 
angular frequency increases. The group velocity of the waves is greater than 
the phase velocity, which shows anomalous dispersion. 
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Surface Plot of Propagation Factor Surface Plot of Attenuation Factor 




Figure 3: Propagation of waves is more near the event horizon. The waves 
grow as they move away from the event horizon. The dispersion is found to 
be anomalous. 
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7 SUMMARY 



This paper is devoted to discussing the isothermal plasma wave properties of 
the Schwarzschild black-hole's magnetosphere by using the 3+1 formalism. 
For this purpose, the GRMHD equations are derived and discussed by taking 
one-dimensional perturbations in perfect MHD flow with its planar analogue. 
These equations are written in component form and then Fourier analyzed by 
using the assumption of plane waves. We have assumed that the black hole's 
gravity does not introduce any effect on magnetic field. The determinant of 
the coefficients are solved to get complex dispersion relations, which yield 
the wave numbers. The properties of plasma are inferred on the basis of 
this number, and the relevant quantities are obtained in graphical form. A 
summary of the results is given below: 

From the three dispersion relations, we find that the wave propagation 
factor becomes infinite. This implies that the wave number is infinite at the 
event horizon; hence, waves vanish there. This corresponds to the well-known 
fact that no information can be extracted from a black hole. 

The attenuation factor assumes extreme values near the event horizon. 
These values gradually decrease or increase as they move away from the 
horizon. The gravity of the black hole causes the attenuation factor to take 
extreme values. Thus, the wave growth or damping is larger near the event 
horizon. The attenuation factor is indefinite at the event horizon. The wave 
propagation factor is high near the black hole's event horizon. As the waves 
go far from the event horizon, the propagation decreases, which shows that 
the perturbations are high near the gravitational well. 

It is worthwhile mentioning here that waves are highly excited near the 
event horizon (Fig. 1). This excitement decreases as they move away from 
the event horizon. This is due to immense gravity of the black hole near the 
event horizon. Similar is the case of Fig. 2, where the waves are randomly 
excited near the event horizon. 

We notice that the dispersion is anomalous here while for a cold plasma 
(Fig. 2) normal wave dispersion was found. This indicates that the 
plasma pressure causes the waves to pass through the region in the neigh- 
borhood of the event horizon. 

The dispersion relations are obtained by using the 3+1 ADM formalism 
and contain the factor of acceleration (depends on the lapse function and is 
equal to —g), which makes them different from the usual MHD dispersion re- 
lations. We have investigated the waves propagating in a plasma influenced 
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by the gravitational field. We observe that the values of k show gravita- 
tional effects on the smooth harmonic wave type perturbations. It would be 
interesting to investigate the plasma wave properties by taking a rotating 
background. 

ACKNOWLEDGMENT 

We appreciate the Higher Education Commission, Islamabad, Pakistan, 
for its financial support during this work through the Indigenous PhD 5000 
Fellowship Program Batch-II. 

References 

[1] J. A. Petterson: Phys. Rev. D10, 3166(1974). 

[2] R. Arnowitt, S. Deser, and C. W. Misner: Gravitation: An Introduction 
to Current Research (Wiley, New York, 1962). 

[3] S. W. Hawking: Nature 248, 30(1974); Commun. Math. Phys. 43, 
199(1975); Phys. Rev. D13, 191(1976). 

[4] J. B. Hartle: Phys. Rev. D8, 1010(1973); ibid. 9, 2749(1974). 

[5] R. S. Hanni and R. Ruflini: Phys. Rev. D8, 3259(1973). 

[6] P. Hajicek: Commun. Math. Phys. 36, 305(1974). 

[7] R.L. Znajek: Ph.D. Thesis (University of Cambridge, 1976); Mon. Not. 
R. Astr. Soc. 185, 833(1978). 

[8] T. Damour: Phys. Rev. D18, 3598(1978); Procs. of Second Marcel 
Grossman Meeting on General Relativity ed. R. Ruffini, North Holland, 
Amsterdam, 587(1982). 

[9] T. Damour: Ph.D. Thesis (University of Paris VI, 1979) 

[10] K. S. Thorne and D. A. Macdonald: Mon. Not. R. Astron. Soc. 198, 
339(1982); ibid 198, 345(1982). 



16 



[11] K. S. Thorne, R. H. Price and D. A. Macdonald: Black Holes: The 
Membrane Paradigm (Yale University Press, New Haven, 1986). 

[12] K. A. Holcomb and T. Tajima: Phys. Rev. D40, 3809(1989). 

[13] K. A. Holcomb: Astrophys. J. 362, 381(1990). 

[14] C. P. Dettmann, N. E. Frankel and V. Kowalenko: Phys. Rev. D48, 
5655(1993). 

[15] R. Khanna: Mon. Not. R. Astron. Soc. 294, 673(1998). 

[16] S. S. Komissarov: Mon. Not. R. Astron. Soc. 336,759(2002). 

[17] X. -H. Zhang: Phys. Rev. D39, 2933(1989). 

[18] X. -H. Zhang: Phys. Rev. D40, 3858(1989). 

[19] T. Regge and J. A. Wheeler: Phys. Rev. 108, 1063(1957). 

[20] F. Zerilli: Phys.Rev. D2, 2141(1970); J. Math. Phys. 11, 2203(1970); 
Phys. Rev. Lett. 24, 737(1970). 

[21] R. H. Price: Phys. Rev. D5, 2419(1972); ibid 5, 2439(1972). 

[22] R. M. Wald: Phys. Rev. D10, 1680(1974). 

[23] J. Sakai and T. Kawata: J. Phys. Soc. Jpn. 49, 747(1980). 

[24] V. Buzzi, K. C. Hines and R. A. Treumann: Phys. Rev. D51,6663(1995); 
ibid 51, 6677(1995). 

[25] K. E. Oughstun and N. A. Cartwright: J. Opt. A: Pure Appl. Opt. 4, 
S125(2002). 

[26] M. Sharif and U. Sheikh: Complex wave numbers in the vicinity of the 
Schwarzschild Event Horizon, Int. J. Mod. Phys. A (2007, to appear). 



17 



